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Abstract 

At high energies, Drell-Yan (DY) dilepton production viewed in the target 
rest frame should be interpreted as bremsstrahlung and can be expressed in terms 
of the same color dipole cross section as DIS. We compute DY cross sections on 
a nucleon target with the realistic parameterization for the dipole cross section 
saturated at large separations. The results are compared to experimental data 
and predictions for RHIC are presented. The transverse momentum distribution 
of the DY process is calculated and energy growth is expected to be steeper 
at large than at small transverse momenta. We also calculate the DY angular 
distribution and investigate deviations from the 1 + cos 2 9 shape. 

PACS: 13.85Qk; 13.85.Lg; 13.60.Hb 

Keywords: Drell-Yan process; dipole cross section; low x 



1 



1 Introduction 



The Drell-Yan (DY) process in the kinematical region where the dilepton mass M is 
small compared to the center of mass energy is of similar theoretical interest as 
deep-inelastic scattering (DIS) at low Bjorken-x. Both processes probe the target at 
high gluon density where one expects to find new physics. In contrast to DIS, where 
only the total cross section can be measured, there is a variety of observables which 
can be measured in the DY process, such as the transverse momentum distribution or 
the angular distribution of the lepton pair. 

The color dipole approach to the DY process suggested by one of the authors U 
(see also 0) provides a convenient alternative to the well known parton model, in 
particular, it is especially appropriate to describe nuclear effects [[[], [|. However, the 
dipole approach was not tested so far for the case of proton-proton collisions. In this 
paper, we calculate different characteristics of the DY process on a proton target. 

Although cross sections are Lorentz invariant, the partonic interpretation of the 
microscopic process depends on the reference frame. In the target rest frame DY dilep- 
ton production should be treated as bremsstrahlung, rather than parton annihilation. 
The space-time picture of the DY process in the target rest frame is illustrated in fig. 
[1]. A quark (or an antiquark) from the projectile hadron radiates a virtual photon on 
impact on the target. The radiation can occur before or after the quark scatters off 



Figure 1: In the target rest 
frame, DY dilepton produc- 
tion looks like bremsstrahlung. 
A quark or an antiquark in- 
side the projectile hadron scat- 
ters off the target color field 
and radiates a massive pho- 
ton, which subsequently de- 
cays into the lepton pair. The 
photon can also be radiated 
before the quark hits the tar- 
get. 

A salient feature of the rest frame picture of DY dilepton production is that at 
high energies and in impact parameter space the DY cross section can be formulated 
in terms of the same dipole cross section as low-x Bj DIS. Note that also the transverse 
momentum distribution of DY dileptons can be expressed in terms of this dipole cross 
section [[| . The crucial input to all calculations is the dipole cross section of interaction 
of a qq pair with a nucleon which at present cannot be reliably calculated. We employ 
the parameterization of Golec-Biernat and Wiisthoff @], since it describes well all DIS 
data in the range of Q 2 which is relevant for DY. 

As a result of the experimental situation, most work in low-x physics is done for 
the case of DIS where the structure function Fi has been measured extensively at 
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HERA. In contrast to this, only very few data are available on the I0W-X2 DY process. 
This will however change with the advent of RHIC. In this paper, we perform the first 
comparison between calculations in the dipole picture of DY and the available data. 
In section |2] we give a description of the color dipole formulation of the DY process. 
The results of the calculations for DY are presented in section |3|. 



2 The DY-process in impact parameter space 

The cross section for radiation of a virtual photon from a quark after scattering on a 
proton, can be written in factorized light-cone form [1], ||, |3|] , 

da(qp^YX) _ f ,2 1 .t.t,l 



din a 



J d 2 p\^ L q (a,p)\ 2 a qij (x2,ap), 



similar to the case of DIS. Here, a qq - is the cross section for scattering a gg-dipole off a 
proton which depends on the qq separation ap, where p is the photon-quark transverse 
separation and a is the fraction of the light-cone momentum of the initial quark taken 
away by the photon. We use the standard notation for the kinematical variables, 
X2 = {■sjxp + 4r — Xp)/2, t = M 2 /s = X1X2, where Xp is the Feynman variable, s is 
the center of mass energy squared of the colliding protons and M is the dilepton mass. 
In ([!]) T stands for transverse and L for longitudinal photons. 

An interesting feature of our approach is the appearance of the dipole cross section 
in ([J), although there is no physical gg-dipole in fig. [|. The physical interpretation 
of (]T|) is similar to the DIS case. The projectile quark is expanded in the interaction 
eigenstates. We keep only the first eigenstate, 

\q) = y/Z~ 2 \q bare ) + ^i L q \ qi *) + ..., (2) 



7*9 

where Z2 is the wavefunction renormalization constant for fermions. In order to produce 
a new state the interaction must resolve between the two Fock states, i.e. they have 
to interact differently. Since only the bare quarks interact in both Fock components 
the difference arises from their relative displacement in transverse plane. If p is the 
transverse separation between the quark and the photon, the 7*g fluctuation has a 
center of gravity in the transverse plane which coincides with the impact parameter 
of the parent quark. The transverse separation between the photon and the center 
of gravity is (1 — a)p and the distance between the quark and the center of gravity 
is correspondingly a p. A displacement in coordinate space corresponds to a phase 
factor in momentum space. The two graphs for bremsstrahlung, where the photon 
is radiated either before or after impact on the target, have the relative phase factor 
— exp(iap • k±), which produces the color screening factor [1 — exp(iap • k±)} in the 
dipole cross section. 

In Born approximation (two gluon exchange) the dipole cross section is indepen- 
dent of energy. The energy dependence is generated by additional radiation of glu- 
ons, which can be resummed in leading ln(l/x) approximation. With help of the 
Weizsacker- Williams approximation and at small separations, the dipole cross section 
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can be expressed in terms of the unintegrated target gluon density, 



47r „ f d 2 k 



1 — exp(ifc^ • p) 



dG{x 2 ,k 2 ± ) 



a qq -(x 2: p) = fa sP 2 j ^ ^ i "dHkiy (3) 

where k± is the transverse momentum exchanged with the target. This is explained in 
some detail in H]. Note that the color screening factor in (|3]) makes the dipole cross 
section vanish like oc p 2 at p — > 0. This salient property of the dipole cross section is 
the heart of the color transparency phenomenon |6|, [7], |8| . 

In terms of Regge phenomenology, the color dipole approach accounts only for 
the pomeron part of the cross section, since the dipole cross section (|3|) is governed 
by gluonic exchange mechanisms. Therefore, this approach can be applied only at 
high energies, i.e. at small x 2 . As already mentioned above, the partonic interpre- 
tation of scattering processes depend on the reference frame. In terms of the parton 
model, which is formulated in the infinite momentum frame of the proton, the dipole 
approach corresponds to annihilation of projectile quarks (antiquarks) with sea anti- 
quarks (quarks) of the target generated via gluons. Note that the statement, whether 
a sea quark belongs to the target or to the projectile, is frame dependent. If the pro- 
jectile quark or antiquark becomes slow in the limit a — > 1, it can be interpreted in 
the infinite momentum frame of the target as anti-seaquark or seaquark of the target 
which annihilates with the projectile parton. No annihilation with valence quarks from 
the target is taken into account in the dipole picture. Note also that valence as well 
as sea parton distributions of the projectile are contained in the parameterization of 
the projectile structure function in ([?]). Therefore, the formulation of the DY process 
presented in this section is not fully symmetric between projectile and target. 

In addition to sea quarks generated from gluon splitting there is a part of the sea 
generated nonperturbatively from the meson cloud of the nucleon ||. This contribution 
has received much attention in connection with the d/u asymmetry measured recently 
by FNAL E866/NuSea JTTJ. Since such a sea component steeply decreases at small x 2 , 
the dipole approach eq. ([l|) can be safely applied in this region. 

The transverse momentum distribution of DY pairs can also be expressed in terms 
of the dipole cross section |J. The differential cross section is given by the Fourier 
integral 

da(qp^^*X) 1 f o , 2 -*mt* / — * *\ t i 

d\nad 2 q ± = (2^ J Pl p2 exp ^ q± ' ~ PV^ri^ PV^T^ PV 



{<j q q{x 2 , api) + (T q g(x 2 , ap 2 ) - a qq -(x 2 , a(pi - p 2 ))} . (4) 



1 

x - 
2 

after integrating this expression over the transverse momentum of the photon, one 
obviously recovers ([!]). The expressions for the LC wavefunctions needed here are 

^ q {a,pi)^ q (a,p 2 ) = ^||mJa 4 Ko(r/pi)Ko(r/p 2 ) 



[1 + (1 - «) 2 ] if^Ki (r) Pl ) Ki ( V p 2 ) }, (5) 
P1P2 J 



^(a,pi)^(a,p 2 ) = -^M 2 (l-a) 2 K (r) Pl )K (r,p 2 ), (6) 



4 



with rj 2 = (1 — a)M 2 — a 2 m 2 . We introduce a quark mass m/ = 200 MeV. The quark 
mass has virtually no influence on the numerical results in pp collisions, fig. |2|, but will 
be more important in proton- nucleus collisions HIT] . Three of the four integrations in 



(||) can be performed analytically for arbitrary a gq - [12 . 

For embedding the partonic cross section ([!]) into the hadronic environment, one 
has to note that the photon carries away the momentum fraction x\ from the projectile 
hadron. The hadronic cross section reads then 



d(7 Gt^iyyx X~\^ 



dM 2 dx F 3ttM 2 x 1 + x 2 



j xi ^22 z f (9/ (-J + 9/ (-) ) — ^ — 



1 f 1 da FP /xi\ dajqp -» 7* A) 
37rM 2 + X2 a 2 v a / d In a 

and similar for the transverse momentum distribution (|j). The factor a em /(37rM 2 ) 
accounts for the decay of the photon into the lepton pair. Remarkably, the parton 
densities qf,qjoi the projectile enters just in the combination .Ff , which is the structure 
function of the proton. Therefore we did not include the fractional quark charge Zf in 
the DY wavefunctions (|5|, ||). The structure function F| is needed at large values of 



XBj- We employ the parameterization from |13[ in our calculations. 

The dipole cross section is largely unknown, only at small distances p it can be 
expressed in terms of the gluon density. However, several parameterizations exist in 
the literature, describing the whole function a q q(x,p), without explicitly taking into 
account the QCD evolution of the gluon density. A very economical parameterization 
is provided by the saturation model of Golec-Biernat and Wiisthoff [Q, 

p 2 QI 



(Tqq(x, p) = (T 



exp 



4(x/x ) x 



where Q = 1 GeV and the three fitted parameters are <r = 23.03 mb, x = 0.0003, 
and A = 0.288. This dipole cross section vanishes oc p 2 at small distances, as implied 
by color transparency and levels off exponentially at large separations, which reminds 
one of eikonalization. The authors of {§] are able to fit all available HERA data with 
a quite low \ 2 an d can furthermore also describe diffractive HERA data. Although 
the parameterization (§) might be unrealistic at very large distances (see discussion in 
TJ|), we can safely use it, because DY data are all taken at quite large virtualities, 



where M] works well. 



3 Calculation of DY cross sections 

We can now proceed and investigate how well DY data are reproduced in the color 
dipole approach. At present, there are however not many data for DY cross sections at 
low X2- We compare to those data for p 2 H scattering from E772 [15] which correspond 
to the lowest values of x 2 . Since the dipole approach is valid at small x 2 , we only 
compare to points with x 2 < 0.1. The result of our calculation, using ([?]), is shown 
in fig. ^. The curves for RHIC (dashed curves) correspond of course to lower values 
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Figure 2: The points represent the measured DY cross section in p 2 H scattering from 
T5fl . Only statistical errors are shown. Note that for these points 0.03 < 22 < 0.09. 
These values of X2 are already quite large for the dipole approach. The curves are 
calculated with the dipole cross section (H) without any further fitting procedure. The 
solid curves are calculated at the same kinematics as the data point (center of mass 
energy y/s = 38.8 GeV). The dashed curves are calculated for RHIC energies, \/s = 500 
GeV. For each energy, the lower curve is for quark mass m/ = 200 MeV, the upper 
curve for rrif = 0. 



for %2 than the ones for E772. The DY cross section increases, because the dipole 
cross section increases with energy. No further fitting procedure of the parameters in 
the dipole cross section @ was performed. The data, and in particular the absolute 
magnitude of the cross section is quite well reproduced, except for few points at low 
mass. We emphasize that the curves in fig. |2| are results of a parameter free calculation. 
Varying the quark mass m/ leaves the numerical results almost unaffected. Note also 
that no i^-factor was introduced. As pointed out above, in pQCD the dipole approach 
corresponds to a resummation of logarithms ln(l/x). However, we do not perform 
a pQCD calculation, but employ the phenomenological parametrization (|8[) which is 
fitted to DIS data. We assume that this parametrization contains also contributions 
beyond the leading-log approximation, as well as nonperturbative effects. Therefore, 
we believe that it is not legitimate to use a i^-factor in our approach. 
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Figure 3: The transverse momentum distribution for DY pairs calculated from at 
xp = 0.625 and M = 6.5 GeV. The curves which flatten at small are calculated 
with the realistic dipole cross section while the other two curves are calculated 
with the small p approximation (H). 

Furthermore, we also calculate the transverse momentum distribution of DY dilep- 
ton pairs from (Q). The result is depicted in fig. |3|. The DY cross section is finite at 
q± = 0, in contrast to the first order pQCD correction to the parton model. In the 
parton model, one has to resum large logarithms log(g_i_/M) from soft gluon radiation 
in order to avoid the divergence at q± = 0. In the dipole approach, the cross section 
does not diverge, because of the saturation of the dipole cross section. In order to find 
out, how sensitive the transverse momentum distribution to the large p-behavior of a q q 
is, we do the same calculation with the small p approximation of (§), 

The result is shown by the dotted curves in fig. |3|. There is no divergence at q± = 0, 
because of the quark mass rrif — 200 MeV, but the cross section at small q±_ is quite 
strongly affected. 

We do not compare to data in fig. Rl because all available data are integrated over 
xf and are therefore contaminated by valence quark contributions. An extraction of 



the low-x part of the transverse momentum distribution is in progress jl6fl . Note that 
the differential cross section increases at large transverse momentum faster with energy 
than at low q±. The reason for this behavior is that at large q± small dipole sizes are 
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Figure 4: The left figure shows the dependence of the parameter A, which describes 
the angular distribution of DY pairs (0) on the dilepton mass. The calculation is 
performed for xp = 0.625. The figure on the right displays the q± dependence of A at 
x F = 0.625 and M = 6.5 GeV. 



predominantly sampled, where the dipole cross section increases more rapidly with 
energy than at large separations. These large separations become more important at 
low q^. The resulting broadening of transverse momenta with energy should not be 
confused with the well known broadening at fixed r = M 2 /s, which can be understood 
from purely dimensional arguments. 

The color dipole approach allows one to calculate separately the cross section for 
longitudinal and transverse photons. Experimentally, different polarizations can be dis- 
tinguished by investigating the angular distribution of DY pairs, which can be written 
as 

— — -ocl + Acos 2 #, 10 

dx F dM 2 dcos6 ' v ; 

where 9 is the angle between muon and the z-axis in the rest frame of the virtual photon. 
The parameter A equals to ±1 for transverse and longitudinal photons respectively. 
Therefore, it can be calculated as 

A = ^^. (11) 

Data for the angular distribution of DY pairs is usually presented in the dilepton center 
of mass frame and the value of A depends on the choice of ^-direction. Since the dipole 
approach is formulated in the target rest frame, it is convenient to put the z-axis in 
direction of the radiated photon ||. The target rest frame and the dilepton center of 
mass frame are then related by a boost in ^-direction. Note that in the dilepton center 
of mass frame, the z-axis is antiparallel to the target momentum. This frame is called 
the w-channel frame and the curves we present for A are valid for this frame. 

We study the dependence of A on the dilepton mass and on the transverse momen- 
tum of the pair. Our results are shown in fig. [|. The deviation of A from unity decreases 
very slowly with increasing mass. Although in fig. [| (left) A is slightly smaller at RHIC 
energies than at E772 energies, the deviation from a 1 + cos 2 9 distribution is typically a 
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5% effect. Note that in DY from pion-tungsten scattering at large xf a, sudden change 
of the angular distribution from 1 + cos 2 9 to sin 2 9 has been observed |17| . This is 



usually explained by interactions of the spectator quark [18] . Such mechanisms are not 
included in the dipole approach. The Kx-part in the transverse light cone wavefunction 
always dominates over the K part in the longitudinal wavefunction. Thus, deviations 
from the 1 + cos 2 9 shape are always small. 

As function of the transverse momentum q±, fig. |] (right), deviations from unity can 
become ~ 10% and A exhibits an interesting nonmonotonous behavior which can be 
checked in future experiments. Note that in the parton model, the Lam- Tung relation 
i~9fl and helicity conservation require that \(q± = 0) = 1, which is obviously not the 



case in the dipole approach. Thus, the Lam- Tung relation is violated in the dipole 
approach. The reason for this behavior is caused by nonperturbative effects, which are 
parameterized in the dipole cross section. Using a parameterization oc p 2 instead of 
(§) would yield X(q± = 0) = 1, as can be seen from eq. (16) in ||. Experimentally, the 
Lam- Tung relation is found to be violated [0, 120 . 

We did not calculate the </> dependence of the cross section, since there is no hope 
that this will be measured within a foreseeable future. The only way to check the 
Lam- Tung relation in the not too far future is to study the limit X(q± — > 0), which 
is possible at RHIC. It is a special virtue of the dipole approach, that one can easily 
perform calculations at q± <C M . 



4 Summary 

In this paper, the first realistic calculations in the color dipole approach to the DY 
process in proton-proton collisions are presented. We employ the parameterization 
[ of the dipole cross section and find good agreement with E772 data [15| at low 



X2, without any i^-factor or free parameter. The quark mass, which is in principle 
undetermined, has virtually no influence on the numerical results. The cross section 
steeply rises with energy and is about four times larger at RHIC than at Fermilab. 

We also study the transverse momentum distribution of DY pairs. As a consequence 
of the saturation of the dipole cross section at large separations, the differential cross 
section does not diverge at zero transverse momentum, in contrast to the first order 
perturbative QCD correction to the parton model. The differential cross section rises 
with energy faster at large than at small momentum transfer. This correlates with 
the fact that the dipole cross section rises with energy steeper at small than at large 
separations. 

We parameterize the angular distribution of the DY pairs as (1 + Acos 2 ^) and 
calculate the coefficient A in the w-channel frame as a function of dilepton mass and 
transverse momentum. We find that A as a function of M is typically around ~ 0.95. In 
addition we find that A does not go to unity for vanishing transverse momentum. This 
is a consequence of the parameterization of the dipole cross section which we employ. 
The behavior of A at small q± can be checked in proton-proton collisions at RHIC and 
possibly could point out the presence of dynamics beyond the QCD improved parton 
model. 
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